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Figure  Captions 


Figure  1 


Naval  Postgraduate  School  R/"/  Acania  at  the  conclusion  of 
CEWCOM-76 . 


Figure  2 


Normalized  heat  conduction  coefficient  (eddy  diffusivity) 
versus  Richardson  number. 


Figure  3 


Variation  of  atmospheric  quantities  at  Z  =  10  meters: 

a)  temperature  structure  fionction,  C 

b)  rate  of  dissipation  of  velocity  fluctuations,  £ 

c)  wind  velocity,  U 


Figure  4 


Figure  5 


Figure  6 


Temperature  structure  function,  C   ,  at  Z  =  10  meters 
versus  time  in  Monterey  Bay.   Solid  circles  are  unstable 
conditions  (12/76)  and   solid  triangles  are  near  neutral 
conditions  (7/76) . 

2 
Diurnal  variation  of  Crp  : 

a)  open  ocean 

b)  coastal  ocean 

R/V  Acania  tracks  off  Southern  California  during  fog  events 
for  CEWCOM-76: 

a)  Sep  27-29     (Fog  #1) 

b)  Oct  5-9     (Fog  #'s  2,3,4) 

The  six  digit  n^miber  indicates  the  day  of  the  month  (dd)  and 
time  of  day  (tttt)  in  the  code  -  ddtttt. 


Figure  7 


Sea  surface  temperature  contours  from  aircraft  IR  measurements 
Courtesy  of  Backes  (1977)  and  Ralph  Markson,  Airborne  Research 
Associates . 


Figure  8 


C   at  Z  =  10  meters  versus  time  for  fog  events 
T 

a)  Fog  #1  Sep  27-28 

b)  Fog  #2  Oct  5 

c)  Fog  #3  Oct  a 

d)  Fog  #4  Oct  9 


Figure  9 


£  at  Z  =  10  meters  versus  time  for  fog  events 

a)  Fog  #1    Sep  27-28 

b)  Fog  #2     Oct   5 

c)  Fog  #3     Oct   8 

d)  Fog  #4    Oct  9 


Figure  10      Sensible  heat  flux,  W,  versus  time  for  fog  events: 

a)  Fog  #1  and  Fog  #2 

b)  Fog  #3  and  Fog  #4 

Figure  11  Radiosonde  measurements  of  temperature  (solid  line) 
and  relative  humidity  (dashed  line)  versus  altitude 
for  fog  events : 

a)  Fog  #1 

b)  Fog  #2 

c)  Fog  #3 

d)  Fog  #4 
Courtesy  of  NWSMET. 


Figure  12 


Ensemble  average  of  turbulence  parameters  one  hour  before 
and  after  fog  is  encountered  and  one  hour  before  and  after 
fog  dissipates: 

a)  Q.     /Q. J-      versus  time 

b)  £/£  versus  time 


2 
Figure  13       Height  dependence  of  C    normalized  to  the  Z  =  6.6  meter 

value  for  fog  (x's)  and  non-fog  (o's)  data.   Neutral 
equilibrium  yields  the  Z~^^  line  indicated,  non-neutral 
equilibrium  would  be  a  less  steep  line  in  this  representa- 
tion . 


Figure  14      Dimensionless  temperature  structure  parameter  versus 

Richardson  number  {atmospheric  stability)  for  fog  (x's) 
and  non-fog  (o's)  data.  The  solid  line  is  the  overland 
results  from  Wyngaard,  et.  al .  (1971). 
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Abstract 

Atmospheric  turbulence  parameters  £  (turbulent  energy  dissipation 

2 
rate) ,  and  C    (temperature  structure  function)  were  measured  in  the 

marine  boundary  layer  during  the  Cooperative  Experiment  for  West  Coast 

Oceanography  and  Meteorology  (CEWCOM-76) .   The  atmospheric  stedsility 

was  characterized  by  measurements  of  profiles  of  mecin  wind  velocity, 

temperature  and  relative  humidity.   This  paper  will  focus  on  turbulence 

properties  and  turbulent  heat  fliox  measured  during  four  meirine  fog 

events.   The  fogs  were  characterized  by  a  slightly  unstable  boundary 

layer  with  moderate  increases  of  heat  flxxx  during  the  fog.  An  ensemble 

2 
average  of  the  fog  events  showed  that  C   emd  £  reached  peak  values 

immediately  before  aind  after  the  fog  was  encountered.   An  cuialysis  of  the 

2 
dimensionless  temperature  structxare  function  (C   normalized  by  the  square 

of  the  temperature  gradient)  gave  lower  values  of  temperature  fluctuations 

in  fog  than  in  clear  air. 


I.   Introduction 

During  September  and  October  of  1976  the  Naval  Postgraduate  School 
(NPS)  shipboard  turbulence  and  profile  systems  were  used  in  a  marine  fog 
and  boundary  layer  research  cruise  off  Southern  California  aboard  the  NPS 
research  vessel  R/V  ACANIA  (Fig.  1).   Termed  the  Cooperative  Experiment  in 
West  Coast  Oceanography  and  Meteorology  (CEWCOM-76) ,  the  cruise  was  jointly 
organized  by  NPS  and  the  Naval  Electronics  Laboratory  Center  (NELC)  of  San 
Diego  in  collaboration  with  CALSPAN  Corporation  under  the  sponsorship  of 
NAVAIR.   Also  represented  on  board  were  the  Naval  Avionics  Facility, 
Indianapolis  (NAFI) ,  the  Naval  Research  Laboratory  (NRL) ,  the  California 
Air  Resources  Board  (CARB) ,  and  a  Naval  Weather  Systems  Service  Mobile 
Environmental  Team  (NWSMET) .   This  was  the  most  recent  of  a  series  of  data 
gathering  cruises  forming  a  major  part  of  a  program  of  marine  fog  investi- 
gation sponsored  by  Naval  Air  Systems  Command.   It  was  also  the  third  such 
cruise  to  include  measurements  of  atmospheric  turbulence.   The  first  was 
conducted  in  1974  aboard  the  R/V  ACANIA  off  Northern  California  (Mack  et  al, 
1975);  the  second  was  conducted  in  1975  aboard  the  USNS  HAYES  off  Novia 
Scotia  (Gathman  and  Larson,  1976). 

Marine  fog  is  one  of  the  least  understood  and  most  complicated  atmos- 
pheric phenomena.   Mack  (1975)  has  identified  at  least  five  different  fog 
formation  "scenarios"  over  the  ocean.   Although  considerable  fog  forecasting 
success  has  been  obtained  by  skilled  phenomenologists  employing  various 
fog  indices  (such  as  inversion  height  and  air  sea  temperature  difference) 
-n  local  coastal  areas  (Leipper,  1948,  Backes ,  1977),  accurate  marine  fog 
prediction  in  the  open  ocean  will  surely  benefit  from  applications  of  know- 
i.euge  from  the  synoptic  to  the  raicroscalfe.   On  the  raicroscale,  the  pro- 
duction of  fog  is  governed  by  several  (interacting)  factors:   aerosols, 


radiation,  hydrostatic  stability,  temperature  inversion  cliaracteristics , 
and  turbulent  transport  of  heat  and  water  vapor.   Turbulent  transport  was 
first  incorporated  into  fog  modeling  by  Rodhe  (1962).   More  modern  fog 
models  (Lee  and  Lin,  1975;  Barker,  1975)  employ  sophisticated  turbulent 
parameter  stability  and  scaling  laws  (see  section  II)  that  are  unverified 
in  the  marine  boundary  layer  (Davidson  et  al,  1977).   The  main  thrust  of 
the  NFS  turbulence  group  during  CEWCOM-76  was  to  attack  these  problems  by: 
1)   providing  modelers  with  measurements  of  turbulence  and  mean  profiles 
in  fog,   2)   obtaining  data  for  evaluation  of  the  turbulence  scaling  laws 
in  the  marine  boundary  layer,   3)   examining  some  of  the  characteristics  of 
turbulence  "background"  levels  for  non-fog  conditions  in  the  marine  boundary 
layer. 
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II.   Theory 

A.   Boundary  layer  equations 

The  boundary  layer  is  that  part  of  the  atmosphere  where  friction 
with  and  heating  by  the  surface  play  an  important  role  in  the  generation 
of  turbulence.   Near  the  surface  the  shear  stress  and  scalar  fluxes  are 
essentially  constant  with  height.   In  this  region  the  fluxes  can  be  rep- 
resented in  terms  of  scaling  parameters  (such  as  u^  and  T^)  that  are  in- 
dependent of  height.   We  shall  refer  to  this  layer  of  nearly  constant  shear 
stress  and  flux  as  the  surface  layer.   For  a  complete  treatment  of  the 
surface  layer  equations  we  suggest  Lumley  and  Panofsky  (1964) ,  Businger 

(1971)  and  Kraus  (1972).   The  normalized  fluxes  of  momentum  (F  ),  heat 

m 

(F,  )  and  water  vapor  (F  )  are 
h  q 

Fjj^  =  -  <'u'w'>  =  u^  (ia) 

F^  =  -  <T'w'>  =  u^  T^  (lb) 

Fq  =  -  Cq'w'>  =  u^  -q*  (Ic) 

where  u'  and  w'  are  the  fluctuating  horizontal  and  vertical  wind  velocities, 
respectively,  T'  is  the  fluctuating  potential  temperature  and  q'  is  the 

fluctuating  specific  humidity.   The  Reynolds  shear  stress,  T,  is  related  to 

2 
the  friction  velocity,  u^,  by  T  =  p  u^  ,  where  P  is  the  density  of  air.  In 

Che  surface  layer  we  can  write  simplified  turbulent  energy  balance  equations 

using  production  terms  (proportional  to  the  vertical  gradient,  3x/3z  j  and 

dissipation  terms  [£  ]  as  follows: 

velocity:   u.   3u/3z  -  ^  u,  T  .  =  £  (2a) 

■'  it  i   *   V* 

temperature:  u^  T^  3T/3z  =  £™  (2b) 

humidity:   u^  q^  9q/3z  =   £  (2c) 
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Note  that  the  velocity  equation  (where  £  =  £)  has  an  extra  pro- 
duction (or  suppression)  term  due  to  the  bouyancy  of  air  (T  is  the  vir- 
tual potential  temperature).  In  the  surface  layer,  the  appropriate  tur- 
bulent transport  scale  parameter  is  the  height  above  the  surface,  z.  In 
this  regime,  turbulence  properties  can  be  scaled  through  the  similarity 
parameter  z/L,  where  L  is  the  Monin-Obukhov  stability  length. 

T  ^* 
L  =  ^^^     ,  (3) 

V* 

where  K  =  .35  is  the  von  Karmon  constant.   We  can  now  define  the  mean  pro- 
file functions  (|)   (z/L)  as  (Businger,et  al,  1971) 

*x  f^/"  =  f  If  (*) 

If  we  define  a  dimensionless  dissipation  function  (Wyngaard  and  Cote,  1971) 
E^  (z/L)  =  — 2  (5) 

then  the  energy  equations  become 

4)^  (z/L)  -  z/L  =  E^  (z/L)  (6a) 

4)^  (z/L)  =  E^  (z/L)  (6b) 

cj>   (z/L)  =  E   (z/L)  (6c) 

q       q 

Present  estimates  of  4*   and  E  are  given  in  Appendix  A. 

In  the  inertial  subrange,   turbulence  is  nearly  isotropic  and  we 
can  represent  the  intensity  of  turbulence  fluctuations  as  a  one-dimensional 
power  spectrum,  S^(k),  of  the  form  proposed  by  Kolmogorov 

S^(k)  =  .25  cj   k"^/^  (7) 
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wherti  k  is  the  wavenumber  and  C    is  Che  structure  function.   If  the 

X 

quantity  x  has  value  x(r)  at  position  r  and  value  x(r  +  d)  at  position 
r  +  d  then 


C   ^   =     ([    x(r)  -  x(r  +  d)]2>  d"^^^  (8) 


2 
In  the  inertial  subrange  C   is  independent  of  d.   Corrsin  (1951)  has  shown 

that  the  structure  function  can  be  related  to  the  dissipation  by 

C  ^  =  4  S  z~^^^  t  (9) 

XX  X  ^  ' 


where  B   is  a  constant,  8^  =  6   =  .81  and  6   =  .52.   We  can  construct  a 
X  T    q  u 

dimensionless  temperature  structure  parameter,  f(z/L),  by  combining  eqs. 
6  and  9  (Wyngaard ,  et  al. ,  1971) 


2     "^  -'^/2 
C^^  =  T^-  z  -^^  f^(z/L)  (10) 


or,  using  eq.  '+ 

C^  =  (aT/3z)2  2^^^  f^(z/L)  (11) 


We  can  also  use  eq.  5  to  obtain  a  similar  expression  for  £. 

£  =  E  (z/L)  (12) 

KZ  ^ 


Since  Che  Monin-Obukhov  length  is  difficult  to  measure,  we  normally  use 
an  alternate  stability  variable-the  gradient  Richardson  number,  R  , 

(3T  /9z) 

^-    =   f-7T-7T-T2  (13) 

1    T   (du/3z) 
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which  is  related  to  z/L  by 


R^  =  (z/L)  (|)^(z/L)/((f)^(z/L))^  (14) 


2 
The  height  dependence  of  C„  and  £  can  be  obtained  from  eqs.  10  and  12;  for 

2    -2/3         -1 
near  neutral  conditions  (R.  ~  0)  C„  -  z     and  e  -    z      .      For  unstable 

2    -4/3 
conditions  (R.  <<  0)  C„  -  z     and  £  -  constant.   For  stable  conditions 

(R.  >  0),  turbulent  eddies  are  increasingly  damped  by  increasing  stable 

stratification  of  the  atmosphere.   Overland  measurements  indicate  complete 

damping  of  turbulence  for  R.  >  .21. 

B.   Flux  calculations 

There  are  four  methods  of  calculating  scalar  fluxes:   eddy  correlation, 

gradient,  turbulence,  and  bulk.   The  eddy  correlation  method  is  a  direct 

measurement  of  x'w'   in  equation  1.   This  requires  a  stable  platform  so 

we  cannot  use  it  for  shipboard  measurements.   The  bulk  method  (Friehe,  1977) 

relates  the   flux  to  the  mean  wind  velocity  and  the  air-sea  temperature  or 

humidity  difference.   We  have  found  this  method  to  be  unsatisfactory  due 

to  difficulties  in  measuring  the  sea  surface  temperature.   The  gradient 

method  relates  the  vertical  gradient  of  the  scalar  to  the  flux  through  the 

eddy  diffusivity,  K  .   For  instance,  the  normalized  heat  flux  is 

F^  =  -  ICj  OT/3z)  (15) 

Unfortunately,  K^  is  a  function  of  stability, 

!C^  =  <  z  u^/0^  (z/L)  (16) 


shown  in  Fig.  2.   Note  that  K^  =  0  for  R  >  .21  in  this  model.   The  tur- 

:^   from  C^' 


2 
bulence  method  is  based  on  calculations  of  T.  from  C„   and  £.   Combining 


eqs.  9  and  6b  we  find 
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T^~  =  K  z  C^'  z^^'^/O.2   u^  .|.^(2/L)) 


(17) 


Substituting  for  u  with  eq.  12,  we  obtain 


F.  = 


1/3 


l£zl 


4/3 


-1/2 


3.2  (E^(z/L))^^^  *^(z/L) 


(18) 


The    sensible  heat  flux  is  related  to  the  normalized  flux  by 


W  =  c  F, 


(19) 


3  3 

where  C  =  1.3x10  joule/ ("C  m  )  is  the  heat  capacity  of  air.   Equations  15 

and  18  have  their  analogous  forms  for  calculation  of  water  vapor  flux.   Un- 
fortunately, instrumental  limitations  and  problems  (see  Section  III)  have 

2 
prevented  us  from  measuring  C   and  Oq/8z)  during  fog  events. 

1 


15" 


III.   InsLrumenCacion  and  AJ^cllysis 

Since  Che  NPS  work  represented  only  a  part  of  the  total  shipboard 
effort  for  CEWCOM-76,  we  will  first  list  the  key   measurements  made  by  the 
other  groups: 

1.  CALSPAN  -  aerosol  spectra  and  chemical  composition,  visibility 
and  profiles. 

2.  NRL  (R.  Jeck)  -  Knollenberg  aerosol  spectrometer. 

3.  NAFI  (J.  Russell)  -  microwave  refractometer. 

4.  CARB  -  ozone,  carbon  monoxide  and  NO  . 

X 

5.  NWSMET  -  twice  daily  radiosondes. 

The  NPS  instrumentation  was  primarily  devoted  to  measurement  of  mean 
and  turbulence  profiles.   A  comprehensive  description  of  this  system  is 
the  subject  of  another  paper  (Houlihan,  et  al.,  1977).   The  profiles  were 
measured  at  four  fully  instrumented  levels  plus  a  sea  surface  temperature 
sensor.   The  level  heights  were  4.2,  6.6,  11.3  and  17.7  meters  above  the 
mean  waterline. 

Mean  temperature  was  measured  with  Hewlett-Packard  Model  HP2801A  quartz 
oscillator  thennometers,  accurate  to  .01  *C.   Mean  humidity  was  measured 
with  Hydrodynamic , Inc .  Digital  II  type  15-1818  LiCl  Dunmore  sensors,  ac- 
curate to  about  3%  RH.   These  humidity  sensors  have  very  slow  response 
above  about  95%  RH  and  response  variations  from  sensor  to  sensor  make  them 
unreliable  under  fog  conditions.   For  this  reason,  we  do  not  have  fog 
event  humidity  profiles  of  sufficient  accuracy  to  calculate  the  humidity 
flux.   The  wind  velocity  was  measured  with  Thomthwaite  Model  101  cup 
anemometers,  accurate  to  about  5%. 

Wind  speed  fluctuations  (for  e)  were  measured  with  Thermo-Systems ,  Inc. 
Model  TSI  1054B  constant  temperature  anemometers  using  TSI  1210  probes  with 

le 


6  rail,  cylindrical  hoc-film  sensors.   The  hot  films  were  calibrated  on  a 

TSI  1125  Calibrator  with  the  cables  actually  used  during  the  shipboard 

2 
measurements.   Temperature  fluctuations  (for  C   )  were  measured  with  GTE 

Sylvania  ac  wheatstone  bridges  (3.0  kHz  carrier)  using  TSI  1210-P.8  micro 

thermal  probes  w:J.th  2.5  u  diameter  platinum  wire.   We  had  planned  to  measure 

2 
water  vapor  fluctuations  (for  C   )  at  one  level  using  an  Electro-Magnetics 

Research  Corporation  Lyman-alpha  sensor,  however  this  instrument  failed 
early  in  the  cruise.   The  mean  data  was  continuously  averaged  and  periodi- 
cally logged  by  a  microprocessor  designated  MIDAS  (Microprogrammable 
Integrated  Data  Acquisition  System)  developed  at  NPS  and  adapted  for  the 
ACANIA  by  Plunkett  (1977).   The  fluctuating  data  was  amplified,  filtered 
and  recorded  on  a  Honeywell  Model  5600  FM  tape  recorder. 

In  addition,  the  NPS  acoustic  radar  (Aeroenvironment  Model  300)  was 
mounted  on  the  ship  and  provided  a  continuous  monitor  of  the  inversion 
height. 

The  analysis  of  the  turbulence  signals  was  auided  by  noise  consider- 
ations.  The  combination  of  long  sensor  cables  (about  150  ft)  and  somewhat 
impure  ship  generated  power  restricted  our  analysis  options — particularly 
for  the  hot  film  systems.   The  velocity  signals  were  processed  with  Fourier 
power  spectrum  analyzers  to  yield  S  (f).   Using  Taylor's  hypothesis 
(k  =  27Tf/u)  and  equations  7  and  9  we  find 


e  =  2.67 


^  [f5/3  s^if)]    ^/2  (20) 


The  spectral  method  has  the  advantage  that  power  source  noise 

2 
(60,120,180  Hz,  etc.)  can  be  ignored.   The  bulk  of  the  C   data  were  ob- 
tained from  paired  sensors  (separated  by  0.3  meters)  using  eq .  8.   In  some 

2 
cases  only  one  sensor  was  operable;  these  were  spectrum  analyzed  and  C 
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was  calculated  using  eq.  7.   Being  based  upon  the  high  frequency  turbulence 
components,  both  of  these  techniques  are  relatively  unaffected  by  ship 
motions . 
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IV.       Results 

A.  General  Background 

Before  plunging  into  the  fog  data,  it  is  of  interest  to  examine 
the  variability  of  turbulence  parameters  in  the  marine  boundary  layer  to 

establish  a  feeling  for  typical  "background"   behavior.   Nominal  variations 

2 

of  C^  and  e  during  a  24  hour  period  (Fig.  3a  and  Fig.  3b)  in  the  open  ocean 

3 
are  relatively  smooth.   Since  e  is  proportional  to  u  ,  we  expect  it  to  be 

2 
more  variable  than  C^  .   A  typical  daily  variation  of  wind  velocity  at 

2 
z  =  10  meters  (u,q)  is  shown  in  Fig.  3c.   Of  course,  C   is  quite  sensitive 

to  changes  in  atmospheric  temperature  gradients  as  we  can  see  in  a  com- 
parison of  measurements  in  Monterey  Bay  under  neutral  conditions  (7/76) 

2 
and  unstable  conditions  (12/76)  in  Fig.  4.   The  diurnal  variation  of  C 

overland  is  dominated  by  solar  heating  of  the  surface  during  the  day  and 

2 
radiative  cooling  at  night.   This  leads  to  a  maximum  of  C   in  the  afternoon 

and  minimum  around  sunrise  and  sunset-   Since  the  ocean  has  both  high  heat 
capacity  and  excellent  thermal  conductivity  (due  to  turbulent  mixing) , 
the  surface  temperature  is  subject  to  only  a  few  tenths  of  degrees  modu- 
lation by  solar  heating  cycles.   In  the  slightly  unstable  (R.  =  -  0.1) 

boundary  layer  we  typically  encountered  in  the  open  ocean,  we  found  that 

2 
radiation  cooling  of  the  atmosphere  produced  a  slight  maximum  of  C   during 

the  night  (Fig.  5a).   In  coastal  areas,  influenced  by  land-sea  breeze  cycles, 

we  found  diurnal  variations  more  characteristic  (but  less  extreme)  of  those 

found  over  land  (Fig.  5b). 

B.  Fog  Data 

This  paper  will  concentrate  on  the  first  four  fog  events  encountered 
during  CEWCOM  76: 

Fog  n  Sept.  27  1100  -  0800  (Sept.  28) 
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Fog  in  Oct.  5  0700  -  1100 

Fog  //3  Oct.  8  1000  -  1800 

Fog  A  Oct.  9  0800  -  1200 

Fog  event  data  are  provided  in  tables  in  Appendix  B.   The  location  of  the 
ACANIA  during  these  events  can  be  established  from  Fig.  6a  (Fog  //I)  and 
Fig.  6b  (Fog  //'s  2,3,  and  A).   Fog  // '  s  1,3,  and  4  were  open  ocean  fogs, 
all  located  west  of  Point  Conception;  Fog  //2  was  a  coastal  fog  located 
just  south  of  Palos  Verdes.   Sea  surface  temperature  contours  (Fig.  7)  for 
the  first  week,  of  October  were  measured  by  Ralph  Markson  in  the  Airborne 
Research  Associate's  aircraft  using  a  Barnes  PRT-5  infrared  radiometer 
(Backes,  1977).   In  order  to  simplify  the  presentation  of  four  levels  of 
turbulence  data,  the  values  from  each  level  have  been  averaged  together 
after  translation  to  a  10  meter  equivalent  value  using  eqs.  10  and  12. 

The  10  meter  equivalent  averages  of  the  temperature  structure  function 

2 
( <C_  >  )  for  the  four  fog  events  are  given  in  Fig.  8.   Note  the  consider- 

^°  2 

able  variability  of  C   compared  to  Fig.  3a  and  Fig.  4.   Figure  9  is  a 

similar  presentation  for  <£>, ^.   In  some  cases  values  of  £  are  not  avail- 
able during  fog  since  the  impact  of  water  droplets  causes  a  noise  spike 
by  cooling  the  film.   In  heavy  fog,  this  effect  totally  masks  the  real 
signal.   We  were  able  to  obtain  £  values  during  Fog  //'s  1  and  4  by  using 
hot  wires  which,  due  to  their  much  smaller  size  (4.5  y  diameter),  are 
less  affected.   The  heat  fluxes  for  the  fog  events  (Fig.s  10a  and  10b) 
were  calculated  from  the  10  meter  averages  using  eqs.  15  and  18.   Rep- 
resentative radiosonde  dew  point  and  temperature  profiles  are  given  in 
Fig.  11. 

The  correlation  of  changes  in  turbulence  parameters  and  fog  events 
can  be  highlighted  by  preparing  "ensemble  averages  of  the  fog  encounters. 
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Eacli  fog  event  was  divided  into  four  one  hour  sections:   one  before  and 
one  after  the  fog  is  first  encountered,  and  one  before  and  one  after  the 
fog  clears   (Fog  //2  was  excluded  from  this  analysis).   If  a  given  fog  event 

was  less  than  two  hours  duration,  the  second  and  third  sections  were  shortened 

2 
to  1/2  the  fog  duration.  Note  that  both  C    and  e  rise  to  a  maximum  about 

30  minutes  before  fog  is  encountered  (Figs.  12a  and  12b).   One  can  also 

2 
examine  the  more  fundamental  behavior  of  C   in  terms  of  fog  and  non-fog 

comparisons  of  height  dependence  (eq.  10)  and  the  dimensionless  temperature 
structure  parameter  (DTSP),  f2(R.),  from  eq .  11.   Although  the  height  de- 
pendence exhibits  no  obvious  disimilarity  (Fig.  13),  the  DTSP  curve  (Fig.  14) 
implies  lower  values  of  C   for  a  given  (9T/3z)  during  fog  events.   This 
is  consistent  with  the  lower  fluxes  obtained  during  fog  by  the  turbulence 
method  as  compared  to  the  gradient  method.   It  is  quite  conceivable  that 
this  is  an  instrumental  effect,  caused  by  the  accumulation  of  liquid  water 
on  the  microthermal  sensors.   Such  an  accumulation,  if  large  enough,  could 
reduce  the  frequency  response  of  the  sensors  sufficiently  to  attenuate  the 
measured  fluctuations.   Based  on  laboratory  measurements  of  frequency  re- 
sponse  in   salt   coated  wires  (Fairall  &  Schacher,  1977)  one  would  expect 

this  effect  Co  be  small  until  the  water  film  was  on  the  order  of  30  microns 

2 
thick.   The  rapid  increase  of  C   before  dissipation  of  fog  and  the  rapid 

decrease  after  dissipation  are  in  the  opposite  direction  of  an  expected 
water  film  effect.   Despite  the  scarcity  of  droplet  impacts  observed  on 
the  130  micron  diameter  hot  films  (implying  the  2.5  micron  diameter  micro- 
thermal  sensors  remained  very  dry)  during  Fog  42,    a  considerable  drop  in 

C^  was  measured. 

i. 

Fog  event  data  from  other  investigators  can  be  found  in  Backes  (1977) 
and  Mack,  et  al .  (1977). 
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V.   Conclusions 

Despite  the  efforts  of  a  number  of  scientists  during  a  month  of  ship- 
board data  gathering,  unequivocal  statements  aibout  the  nature  of  marine 
fog  are  still  difficult  to  come  by.   Although  it  might  be  wise  to  resist 
the  temptation  to  embark  into  arm  waving  arguments  about  the  reasons  for 
fog  formation,  it  is  interesting  to  note  that  the  three  open  ocean  fogs 
described  "were  found  in  the  vicinity  of  the  cold  water  intrusion  at  Point 
Conception  (upper  left  hand  corner  of  Fig.  6) .   To  describe  these  fogs 
as  being  formed  by  "warm  air  blowing  over  cold  water"  is  probably  much 
too  optimistic,  since  all  were  characterized  by  unstable  surface  layers. 
Indeed,  the  vertical  heat  flxixes  often  dramatically  increased  during  the 
fog.   We  can  make  some  statements  about  the  nature  of  turbulence  in  fog, 
keeping  in  mind  that  they  aire  subject  to  questions  about  the  influence 
of  water  droplets  on  sensors .   The  temperature  structure  function  (and 
therefore  the  turbulent  heat  transport)  was  subject  to  considerably  more 

variations  in  the  atmospheric  conditions  conducive  to  fog  formation  thaxi 

2 
in  the  typical  non-fog  background.   An  ensemble  analysis  of  C   and  £ 

showed  increased  tiarbulence  levels  just  before  and  just  after  the  fog  was 

encoxintered .   This  suggests  that  the  fog  is  often  superimposed  on  some 

large  scale  cellular  structure  in  the  boundary  layer  (on  the  order  of  10 

km  in  the  down  wind  direction) .   When  analyzed  in  dimensionless  form 

2 
(normalized  by  temperature  gradient) ,  C   is  lower  in  fog  than  out  of  fog. 

It  is  possible  the  fog  droplets  are  retarding  the  temperature  fluctuations 

by  evaporation-condensation  response . 

In  general ,  we  have  made  considerable  progress  towards  the  specific 

goals  enumerated  in  the  introduction.   Notably  in  the  provision  of  a 
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body  of  data  for  modelers  to  work  with.  However,  there  is  considerable 
room  for  improvement  (particularly  in  the  water  vapor  flux  capacity)  in 
future  work. 
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Appendix  A 

The  forms  of  the  mean  profile  functions  are  taken  from  Businger,  et  al  (1971) 
.|)^(z/L)  =  (1  -  15  z/L)"^^^         z/L  <  o 

4)^(z/L)  =  (1  +  4.7  z/L)  z/L  >  o 


)^(z/L)  =  .74  (1-9  z/L)"-^''^      z/L  <  o 


0^(z/L)  =  (.74  +  4.7  z/L)  z/L  >  o 


The  dimensionless  velocity  dissipation  function  is  from  Wyngaard  and  Cote 
(1971). 

E^(z/L)  =  [1  +  .5|z/l|-^^]^^^      z/L  <  o 
E^(z/L)  =  [1  +  2.5(z/L)^^^]-^^^     z/L  >  o 

The  dimensionless  temperature  structure  functions  is  from  Wyngaard,  et  al 
(1971). 

f^(z/L)  =  4.9  (1-7  z/L)"^^^      z/L  <  o 

f^(z/L)  =  4.9  (1  +  2.75  z/L)       z/L  >  o 
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